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Indirect CP violation in the D° system can be probed by measuring the parameter ,4r, defined as 
the CP asymmetry of the effective lifetime of the D° meson decaying to a CP eigenstate. This can 
be significantly enhanced beyond Standard Model predictions by new physics. Measurements 
of Af using D° —> K K and D°—>• tt 1 n decays reconstructed from pp collisions collected by 
the LHCb experiment, corresponding to an integrated luminosity of 1.0 fb 1 , are presented. The 
results are 


A r (7T7t) = (+0.33±1.06±0.14) x 10~ 3 , 

A r (KK) = (—0.35 ±0.62±0.12) x 10^ 3 , 

where the uncertainties are statistical and systematic, respectively. These are the most precise 
measurements of their kind to date, and show no evidence of CP violation. 
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Figure 1: Fits to (left) the D° invariant mass distribution and (right) the Aw; = m(D ' 1 ) — w(D°) distribution 
for D° -> K + K candidates from the data subset with magnet polarity down, recorded in the earlier of the 
two running periods. 

1. Introduction 

Similarly to the B° and B* 1 systems, the mass eigenstates of the D° system, IDj^), with masses 
mp 2 and widths r^, are superpositions of the flavour eigenstates |Di, 2 ) = p|D°) ± g|D°), where p 
and q are complex and satisfy \p\ 2 + \q\ 2 = 1. This causes mixing between the |D°) and |D°) states, 
and allows for “indirect” CP violation in mixing, and in interference between mixing and decay, 
when decaying to a CP eigenstate. Indirect CP asymmetries in the D° system can be significantly 
enhanced beyond Standard Model (SM) predictions by new physics [1], In decays of D° mesons 
to a CP eigenstate /, indirect CP violation can be probed using [2] 

^(A m +A d )ycos(j)-xsiruj) , 

where f is the inverse of the effective lifetime of the decay, r/cp is the CP eigenvalue of /, 

x = 2 (m 2 -mi)/(ri +r 2 ), y= (r 2 -r 1 )/(r 1 +r 2 ), a,„ = ( \q/p\ 2 -\p/q\ 2 )/(\q/p \ 2 + \p/q\ 2 ), 

Ad = ( \Af\ 2 — |Ay| 2 )/(|Ay| 2 + |A/| 2 ), with A f the decay amplitude, and (j) = arg(qAf/pAf). The 
effective lifetime is defined as the average decay time of a particle with an initial state of |D°) or 
|D°), i.e. that obtained by fitting the decay-time distribution of signal with a single exponential. 

The LHCb detector at the LHC, CERN, is a forward-arm spectrometer, specifically designed 
for high precision measurements of decays of b and c hadrons [3], During 2011 the experiment 
collected pp collisions at y/s = 7 TeV corresponding to an integrated luminosity of 1.0 fb . Due 
to the large cc production cross section [4], the decay-time resolution of approximately 50 fs for 
D° decays [5] and the excellent separation of n and K achieved by the detector [6], it is very well 
suited to measure Ap with high precision. 

2. Methodology 

The decay chain D* + —> D°7t s f is used to determine the flavour of the D° candidates at pro¬ 
duction, via the charge of the n s meson. The CP-even K K~ and 7t 1 7i ~ final states are used to 
calculate A\ [7]. The predominant candidate selection criteria require the K + K” or 7t + 7t tracks to 
have large impact parameter (IP), large transverse momentum (pp), invariant mass within 50 MeV 


Ar = 


r(D u -> /) - r(D u -> /) 

f (D° —> /) + f (D° -A /) 


V cp 
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Figure 2: Fits to (left) the D° decay-time distribution and (right) the ln(xf p ) distribution for D° —> K + K 
candidates from the data subset with magnet polarity down, recorded in the earlier of the two running periods. 

of the world average D° mass, and for the vector sum of their momenta to point closely back to 
the position of the pp collision. Using data corresponding to an integrated luminosity of 1.0 fb 1 , 
4.8M D °—> K K candidates and 1.5M D °—> n + n~ candidates are selected. The data are divided 
by D° flavour, the polarity of the LHCb dipole magnet, and two separate running periods. Combi¬ 
natorial and partially reconstructed backgrounds are discriminated using a simultaneous fit to the 
distributions of D° mass and Am = m(D* + ) — ra(D°). Examples of these fits are shown in Fig. 1 
for D° —> K K candidates, for data recorded with the magnet polarity down during the earlier of 
the two running periods. 

A tit to the decay-time distribution of the candidates is then used to determine the effective 
lifetimes of the D° and D° signal. Only candidates for which the D* + is produced directly at the 
pp collision are considered as signal. The background from B —y D* X decays is discriminated by 
simultaneously fitting the distributions of the decay time and the natural logarithm of the % 2 of the 
hypothesis that the D° candidate originates directly from the pp collision (ln(;£ IP )). The selection 
efficiency as a function of decay time is obtained from data using per-candidate acceptance func¬ 
tions, as described in detail in Ref. [8]. The decay-time and ln(^f p ) distributions for combinatorial 
and specific backgrounds are obtained from the data using the discrimination provided by the mass 
and Am fits to employ the s Weights technique [9] with kernel density estimation [10]. Figure 2 
shows fits to the distributions of decay time and ln(^ T 2 p ) for D° — > K 1 K candidates, using the 
same data subset as Fig. 1. Inaccuracies in the fit model are examined as a source of systematic 
uncertainty, as discussed in the following section. 

3. Results and systematics 

The fits detailed in the previous section find 

A r (7T7t) = (+0.33±1.06±0.14) x 10~ 3 , 

A r (KK) = (—0.35 ±0.62 ±0.12) x 10 3 , 

where the uncertainties are statistical and systematic, respectively. These are the most precise 
measurements of their kind to date, and show no evidence of CP violation. The dominant systematic 
uncertainties arise from the modelling of the selection efficiency as a function of decay time, and 
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Figure 3: The world averages of (left) Ap and (right) direct vs. indirect CP violation in D° —> h + h decays, 
reproduced from [11]. 


the modelling of the background from B —> D* + X decays. Figure 3 (left) shows the world average 
of A|% which is dominated by these measurements and is consistent with zero. Figure 3 (right) 
shows the combined fit to measurements of direct and indirect CP violation in D° —> h h decays, 
which yields a p-value for zero CP violation of 5.1 % [11]. 

The precision of these measurements will be improved by the addition of 2.1 fb 1 of data 
already collected during 2012. Together with data to be recorded in run II, and, in time, following 
the LHCb upgrade, measurements with precisions of approximately 1 x 10 4 are possible, giving 
great potential for the discovery of indirect CP violation in the D° system. 
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